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ABSTRACT 


This  report  describes  the  progress  accomplished  in  the  past 
quarter  on  the  direct-scattering  computer  program,  its  anticipated 
capabilities  and  limitations,  the  mathematical  description  of  the  target, 
the  physical-optics  formulation,  and  recommendations  for  subsequent 
computer  program  improvement  effort. 
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I.  GOALS 


Our  purpose  is  to  develop  new  theoretical  and  computational  tech¬ 
niques  for  electromagnetic  scattering.  The  end  result  will  be  a  digital 
computer  program  designed  to  calculate  the  scattering  properties  of 
a  wide  class  of  targets. 


II.  INTRODUCTION 

Obviously  there  would  be  many  applications  for  a  computer  program 
that  would  permit  rapid  and  accurate  calculations  of  the  scattering  pro¬ 
perties  of  every  conceivable  target  size,  shape  and  material  in  the  bistatic 
situation  with  an  incident  plane  wave  having  any  given  polarization, 
propagation  direction ,  and  waveform.  However,  the  current  state  of  the 
art  in  scattering  theory  and  computer  technology  is  such  that  high  accuracy, 
speed  and  generality  are  mutually  exclusive.  Therefore,  our  objective 
is  to  develop  a  computer  program  with  suitable  tradeoffs  to  permit  rea¬ 
sonably  accurate  calculations  with  as  much  speed  and  generality  as  is 
feasible.  Since  the  computer  program  will  be  applied  primarily  to  problems 
where  high-frequency  scattering  theories  are  appropriate,  such  techniques 
as  physical  optics  and  the  geometrical  theory  of  diffraction  are  most 
suitable. 

The  physical  optics  technique  yields  reasonably  accurate  scattering 
data  in  many  cases,  and  yet  it  permits  rapid  calculations  for  a  wide 
class  of  targets.  Furthermore,  improved  accuracy  can  be  obtained 
by  making  a  fairly  simple  modification  to  the  physical  optics  solution 
(to  remove  the  erroneous  contribution  from  the  shadow  boundary)  and 
including  the  effects  of  creeping  waves  which  travel  around  the  shadowed 
region  of  the  target. 

The  following  sections  in  this  report,  describe  our  progress  with 
the  creeping-wave  approach  during  the  past  quarterly  period,  set  forth 
the  probable  capabilities  and  limitations  of  the  final  computer  program, 
and  present  the  equations  for  the  physical-optics  formulation.  Finally, 
recommendations  are  given  for  subsequent  improvements  that  should 
be  made  in  the  computer  program  toward  greater  accuracy  and  wider 
applicability. 


III.  PROGRESS  WITH  THE  CREEPING-WAVE  APPROACH 

A  computer  subroutine  for  wedge  diffraction  coefficients  was 
obtained  from  another  group  at  this  Laboratory.  This  program  com¬ 
putes  the  far-zone  scattered  field  of  a  wedge  with  arbitrary  included 
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angle.  The  logic  for  determining  which  wedges  on  the  target  contribute 
to  the  back-scattered  field  was  written  but  has  not  yet  been  tested. 

The  subroutines  for  determining  the  radii  of  curvature,  differential 
arc  length,  and  Gaussian  curvature  at  an  arbitrary  point  on  the  target 
have  been  written  and  tested  for  spheres,  spheroids  and  ogives. 

The  numerical  integration  subroutine  for  the  path  lengths  of  the 
creeping  waves  has  been  written.  The  subroutine  for  determining  the 
specular  point  and  the  attachment  and  reradiation  points  of  the  creeping 
wave  is  being  tested. 

These  subroutines  will  next  be  combined  to  compute  the  backscatter 
from  bodies  of  revolution  for  parallel  polarization.  The  case  of  per¬ 
pendicular  polarization  will  require  further  development  of  a  subroutine 
for  determining  the  creeping  wave  paths. 

The  theoretical  basis  for  the  creeping-wave  approach  is  presented 
in  Reference  1. 

IV.  COMPUTER  PROGRAM  CAPABILITIES 

This  section  sets  forth  the  anticipated  capabilities  and  limitations 
of  the  computer  program  which  will  be  the  end  product  of  this  effort. 

Some  of  the  limitations  could  be  eliminated  with  the  one -year  follow¬ 
up  program  recommended  in  Section  VII. 

Figure  1  shows  a  general  flow  chart  for  the  computer  program. 

The  program  will  handle  perfectly  conducting  targets.  The  surface 
is  considered  to  be  convex  and  continuous.  The  target  may  have  any  size, 
but  the  accuracy  of  the  output  data  will  be  doubtful  when  the  target  is  too 
small  in  comparison  with  the  wavelength,  and  the  computation  time 
and  memory  storage  requirements  will  increase  for  large  targets. 

The  effects  of  antennas  mounted  on  or  in  the  target  will  not  be 
calculated.  The  effects  of  fins  will  not  be  considered,  with  the  possible 
exception  of  fins  on  a  body  of  revolution  for  the  CW  backscatter  case. 

Bistatic  scattering  data  will  be  calculated  with  the  physical -optics 
technique  for  the  CW  and  pulse  cases.  In  the  pulse  case,  there  will 
be  no  limitations  on  the  pulse  length  aside  from  the  requirement  that 
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Fig.  1.  General  flow  chart  for  the  computer  program. 
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the  incident  pulse  have  M  complete  cycles  (where  M  is  an  integer). 

The  incident  wave  is  considered  to  be  an  amplitude -modulated  carrier 
with  a  square-wave  modulation  envelope. 

The  creeping  wave  contributions  will  be  calculated  to  yield  greater 
accuracy  for  CW  backscatter  from  bodies  of  revolution.  The  off- 
diagonal  elements  in  the  scattering  matrix  will  be  set  equal  to  zero  in 
backscatter  problems,  since  a  better  solution  is  not  feasible  at  this 
time. 


The  direct-scattering  computer  program  will  not  calculate  the 
Doppler  frequency  shifts  associated  with  target  motion,  coordinate 
transformations  between  the  target-fixed  reference  frame  and  other 
reference  frames,  or  the  time  and  phase  delays  associated  with  the 
distances  from  the  target  to  the  transmitting  and  receiving  antennas. 


V.  TARGET  DESCRIPTION 

For  physical-optics  calculations  ,  the  target  will  be  described  in 
terms  of  a  large  number  of  points  (xj,yi,z^)  on  the  surface.  These 
points  will  lie  on  the  intersections  of  the  target  surface  with  the  planes 
z  =  Zj  ,  z  =  z2  ,  etc.  The  coordinate  origin  will  be  located  in  the 
interior  of  the  target.  The  target  surface  will  be  approximated  by 
triangular  facets  with  vertices  at  the  given  points.  It  will  be  assumed 
that  the  points  will  cover  the  surface  with  a  density  such  that  each  facet 
is  small  in  comparison  with  the  wavelength. 

The  creeping-wave  analysis  will  be  programmed  for  backscatter 
data  for  perfectly  conducting  bodies  having  axial  symmetry  with  respect 
to  the  z  axis.  The  coordinate  origin  will  be  located  in  the  interior  of 
the  target.  The  surface  will  have  N  analytic  subsurfaces  with  boundaries 
defined  by  the  intersections  with  the  cones  0  =  Subsurface  I  will  be 

described  by  a  set  of  six  real  coefficients  A1(I)S  A2(I)  ,  ...A6(I)  in 
accordance  with  the  following  equation: 


( 1)  F(p  ,z)  =  A  1(1)  p  2  +  A2(I)  z2  +  A 3(1)  p  z  +  A4(I)  z  +  A5(I)  p  +  A6(I)  =  0 

where  p  denotes  the  radial  coordinate  in  the  cylindrical  system.  Thus, 
the  trace  of  each  subsurface  on  the  xz  plane  will  be  a  conic  section. 
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VI.  THE  PHYSICAL  OPTICS  FORMULATION 


If  the  transmitting  antenna  is  at  a  great  distance  from  the  target, 
it  will  illuminate  the  target  with  an  incident  field  which  is  essentially 
a  plane  wave.  In  the  CW  case  we  let  the  time  dependence  eJwt  be 
understood  and  represent  the  incident  electric  field  intensity  as  follows 


where  ( r^ ,  »  4*i)  are  the  spherical  coordinates  of  the  transmitting 

antenna,  (r^,  4i)  are  the  corresponding  unit  vectors, 


(3)  k  =  2-rr/X. , 


\  denotes  the  wavelength,  and  Eg  and  E^  are  complex  constants.  An 
arbitrary  point  on  the  target^si^’face  is  assigned  the  coordinates 
(r,0,4>),  the  unit  vectors  (r,  0 ,  <j> ) ,  and  the  position  vector 

A 

(4)  r  =  r  r 


Finally,  the  receivin^an^en^a  is  assigned  the  coordinates  (rs,0s,<t>s) 
and  the  unit  vectors  (  rs  >  es  » 4>s)  •  Thus  &  and  4>i  specify  the  incidence 
angles  and  0g  and  4>s  are  the  scattering  angles.  The  magnetic  field 
intensity  of  the  incident  plane  wave  is  given  by 


(5)  h‘  =  (tj  Ej,  4i  e‘>  ! — 

T1 


wh  ere 


(6)  ri  =  N/p/e 


The  current  density  induced  on  the  illuminated  portion  of  the  target 
surface  is  approximated  as  follows: 
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(7)  J  =  2$i  X  H1 


where  ft  denotes  the  outward  unit  normal  vector  on  the  surface.  The 
vector  potential  for  the  distant  scattered  field  is  given  by 


(8) 


A  = 


H- _ 

4lT  rs 


ds 


At  a  great  distance  from  the  target,  the  scattered  field  is 


(9)  E  =  -  jcoA  = 


JU)^  -jkr. 


2  it  r . 


i  Jk $  •  r 
e  s  —  ds 


From  Eqs.  (5)  and  (9), 


10)  E 


s  e“jkrs  a 


<6i  EJ>  '  A  EQ>  x  s 


where 


11)  S  =  (jA)  \  ^  ft  ejk(ri  +  rs) '  -  ds 


The  distant  scattered  field  is  represented  by 


(12)  Es  =  $s  Ee  +  $s  Ej) 


-jkr, 


s.  e  J  s 
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s  s 

where  E  and  E,  denote  complex  constants.  From  Eqs.  (10)  and  (12) 
and  the  following^  vector  identity, 


(13)  A  •  (B  X  C)  =  (A  X  B)  •  C 


it  is  found  that 


04)  Eg  =  (Ei  X  -&s  +  E^fc.X  0i}  ■  S 
and 

(15)  e£  =  (E;  +  E;  %bx\)  *  S 

It  is  convenient  to  define  the  CW  scattering  matrix  as  follows: 

(E 

(16) 

lE 

From  Eqs.  (14)  through  (16),  the  i  omplex  elements  in  the  scattering 
matrix  are  given  by 

(17)  Sn  =  (4>i  X  fcs)  •  S 

(18)  S12  =  X  \)  ■  S 

(19)  S21  =  ($t  X  %8)  ■  S 

(20)  S22  =  (4>s  X  0.)  •  S 
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It  is  convenient  to  define  a  "vector  area  function"  A(w)  as  follows 


(21) 


S  =  (j/M  ^  Al 


A(w) 


jpw 


dw 


where 


(22) 


w  = 


/A  ,  A  , 
r.  +  r  )  •  r 
1  s  — 


A  ,  A 
r-  +  r 
i  s 


and 


(23)  p  =  k 


A 

r 


It  may  be  seen  from  Eq.  (22)  that  w  represents  one  of  the  coordinates 
of  a  point  on  the  target  surface,  in  a  rectangular  coordinate  system 
that  is  rotated  in  space  with  respect  to  the  (x,y,z)  system.  The  w  axis 
is  coplanar  with  r •  and  rs  and  bisects  the  angle  between  these  unit 
vectors . 

For  any  point  (x,y,z)  on  the  target  surface, 


(24)  w 


r 

x(sin0-  coscj).  +  sin9  coscf>  )  +  y(sin0-  sincf).  +  sin0  sincf)  ) 
11  SS  11  ss 


+  z(cosO.  +  COS0  ) 

'  1  s' 
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lA 

lri 


where 


(sin0,  cos<{>.  +  sinO  cos4>  )2 
11  s  ® 


+  (sin0.  sin<|>.  +  sin0  sin6  )2  +  (cos0.  +  cos©  )2 
1  1  s  s  1  s 


Once  the  vector  area  function  has  been  calculated,  S  can  be  determined 
efficiently  from  Eq.  (21)  by  numerical  integration.  Equation  11  would 
take  more  computation  time  since  it  involves  a  surface  integral  instead 
of  a  line  integral. 

Once  S  has  been  calculated,  the  elements  in  the  scattering  matrix 
are  determined  as  follows: 


(27)  S12  =SxX12  +Sy  Y12  +Sz  Z12 


(28)  S21  =  Sz  Z21 


(29)  S22  =Sx  X22  +Sy  Y22 


where 


(30)  S  =  x  Sx  +  y  Sy  +  z  Sz 


(31)  Xn  =  -  cos^  sin0g 


(32)  Yu  =  -  sin4>.  sin0g 


(33)  Zu  =  -  (sin(j)i  sin<t>s  +  cos4>.  cos<t>s)  cos0g 


(34)  x12  =  cos0i  sin«t>i  sin0g  -  sin0i  cos0g  sin4>g 
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(35)  Y12  =  sin0i  cos0s  cos$s  -  cos0i  cos^  sin0g 

(36)  Z12  =  (sin4>i  cos4>g  -  cos4>.  sin4>s)  cos@i  cos0g 

(37)  Z2X  =  c.os4>i  sin4>s  -  sin^  cos<t>s 
(  38)  X 22  =  -  sin0i  cos<j)s 

( 39)  Y22  =  -  sin0i  sin<j>g 

(40)  Z22  =  -  (cos^  cos<t>s  +  sin<j>i  sin<t>g)  cos0i 

In  the  pulse  case,  the  incident  wave  is  considered  to  have  M 
complete  cycles  and  a  square  modulation  envelope.  M  is  assumed  to 
be  an  integer,  and  the  incident  field  is  expressed  by 

(41)  E1  =  Eg  +  %i  E*)  sin  (cot  +  k  ri  •  r) 

u(cot  +  k  ri  •  r)  -  u(cot  +  k  ri  .  r  -  cot) 


where  Eg  and  E^j,  are  real  constants,  u(x)  denotes  the  unit  step  function, 
f  is  the  carrier  frequency, 


(42)  co  =  2  it f  , 

(43)  T  =  1/f, 

and 

(44)  t  =  MT  . 
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The  scattered  field  is  given  by 


(45)  ES  (rs,t)  =  [  0S  E®(t  -  rs/c)  +  $g  E®  (t  -  ra/c)]  j- 


It  is  convenient  to  write  the  following  matrix  equation 


(46) 


/f„(  t) 

F12  (t)^ 

(Ee\ 

\f2i  (t) 

fSJ 

fM 

Ua 

The  pulse  response  of  the  target  is  thus  defined  with  four  func- 
tions  of  time  given  by: 

(47)  F„(t)  =  ($.  X  eg)  .  F  (t) 

(48)  F j2  (t)  =  (0g  X  0.)  •  F  (t) 

(49)  F21  (t)  =  X  $g)  •  F  (t) 

(50)  F 22  (t)  =  ($g  X  0.)  •  F  (t) 
where 


(51)  F  (t)  =  £  j  1  ncos  M  +  (3 w) 

‘  [u(wt  +  pw)  -  u(wt  +  P w  -  ut)]  ds 

Eq.  (51)  can  also  be  written  as  follows 


1  1 


(52)  F(t)  =  -  A  (w)  cos  (cot  +  (3 w)  [u(cot  +  pw)  -u(u>t  +  Pw  -  u  t)  ]  dw 

Equation  (52)  will  be  employed  in  the  computer  program  since  it 
permits  more  rapid  calculations  than  Eq.  (51). 

The  scattering  waveforms  Fn(t),  F12(t),  F21(t)  and  F22(t)  will 
be  calculated  at  a  discrete  set  of  equally  spaced  points  in  time. 

The  spacing  between  these  points  will  be  determined  by  an  integer  L, 
supplied  by  the  operator,  as  follows 


(53)  At  =  T/L 


The  first  and  last  points  will  coincide  with  the  initiation  and  the 
termination  of  the  pulse  response. 


VII.  RECOMMENDATIONS 

A  one-year  follow-up  effort  is  recommended  to  increase  the 

accuracy  and  widen  the  scope  of  the  direct-scattering  computer 

program.  Our  recommendations  are  listed  below. 

1.  The  program  should  be  extended  to  handle  imperfectly  conducting, 
absorbe r -coated ,  dielectric  and  composite  targets.  This  could 
be  accomplished  with  the  formulation  developed  in  Reference  2. 

2.  The  scattering  effects  of  antennas  mounted  on  the  target  should 
be  incorporated  in  the  program.  This  is  believed  to  be  feasible 
if  the  input  data  include  the  antenna  input  impedance,  the 
antenna  pattern,  and  the  load  impedance. 

3.  The  capability  of  handling  fins  and  concave  regions  on  the  target 
should  be  developed  by  using  multiple -bounce  reflection  theory. 

4.  The  creeping-wave  subroutines  should  be  generalized  to  the 
bistatic  case  and  to  targets  having  no  axis  of  symmetry. 

5.  An  attempt  should  be  made  to  program  an  accurate  impulse 
response  for  targets  of  arbitrary  shape.  This  might  be  accom¬ 
plished  by  starting  with  the  physical-optics  step  response 
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(calculated  in  the  present  program)  and  applying  various  correc¬ 
tions  with  the  techniques  developed  by  Kennaugh  and  Moffatt.  3 
This  impulse  response  could  then  form  the  basis  for  accurate 
scattering  data  even  when  the  incident  radar  pulse  has  frequency 
modulation  (as  with  chirp  radar)  or  amplitude  modulation  with 
an  envelope  more  realistic  than  the  square  wave. 
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